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We investigate biases in cluster ages and [Fe/H] estimated 
from the (V — K) — (V - I) diagram, arising from inconsis- 
tent Horizontal Branch morphology, metal mixture, treatment 
of core convection between observed clusters and the theoret- 
ical colour grid employed for age and metallicity determina- 
tions. We also study the role played by statistical fluctuations 
of the observed colours, caused by the low total mass of typ- 
ical globulars. Synthetic samples of globular cluster systems 
are created, by means of Monte-Carlo techniques. Each sam- 
ple accounts for a different possible source of bias, among the 
ones addressed in this investigation. Cumulative age and [Fe/H] 
distributions are then retrieved by comparisons with a refer- 
ence theoretical colour-colour grid, and analyzed. Horizontal 
Branch morphology is potentially the largest source of uncer- 
tainty. A single-age system harbouring a large fraction of clus- 
, ters with an HB morphology systematically bluer than the one 
{SJ ' accounted for in the theoretical colour grid, can simulate a bi- 
, modal population with an age difference as large as ~8 Gyr. 
Jj^ ' When only the redder clusters are considered, this uncertainty 
, is almost negligible, unless there is an extreme mass loss along 
the Red Giant Branch phase. The metal mixture affects mainly 
the redder clusters; the effect of colour fluctuations becomes 
negligible for the redder clusters, or when the integrated My 
is brighter than ~ -8.5 mag. The treatment of core convection 
is relevant for ages below ~4 Gyr. The retrieved cumulative 
, [Fe/H] distributions are overall only mildly affected. Colour 
fluctuations and convective core extension have the largest ef- 
fect. When lcr photometric errors reach 0.10 mag, all biases 
found in our analysis are erased, and bimodal age populations 
with age differences of up to ~8 Gyr go undetected. The use of 
both (U -I)-(V - K) and (V - I) - (V - K) diagrams may 
help disclosing the presence of blue HB stars unaccounted for 
Q ■ in the theoretical colour calibration. 

^ ■ Key words. Galaxies: evolution - Galaxies: formation 
. . | Galaxies: star clusters - globular clusters: general 
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Abstract. Age and metallicity estimates for extragalactic globular clusters, from integrated colour-colour diagrams, are exam- 
ined. 



1. Introduction 

A large number of studies about extragalactic globular cluster 
(GC) systems has appeared in the last decades (see, e.g. the re- 
views by Ashman & Zepf fTWSl and Brodie & S trader I2006V 
The main motivation for these analyses is that GCs are ex- 
pected to be good tracers of major episodes of star formation 
of spheroids (see, e.g., Schweitzer l 19971 Kissler-Patig, Forbes 
& Minniti fl998t and in fact inferred GC age/metallicity distri- 
butions have been widely used to investigate the formation and 
evolution histories of early-type galaxies (see, e.g., Kuntschner 
et al. l20Q2l Larsen et al. 120011 Strader & Brodie 120041 Yi et 
al. 2004 Hempel et al. 2005 and references therein). One of 
the most relevant discoveries in the field is the fact that most (if 
not all) extragalactic GC systems show a (integrated) colour bi- 
modality (see, e.g., Gebhardt & Kissler-Patig 1999 1 indicating 
two distinct GC subpopulations, that is, two major star-forming 
episodes or at least two different formation mechanisms in the 
galaxy evolution histories (see Yoon et al. 2006 for an alterna- 
tive view on this issue). 

The interpretation of this bimodality in terms of age 
and/or chemical composition differences is complicated by the 
well known age-metallicity degeneracy (e.g. Worthey 1994) 
whereby the effect on the integrated colours of an age increase 
at fixed metallicity, can be reproduced also by suitably in- 
creasing the metallicity while keeping the age fixed. The age- 
metallicity degeneracy can be broken in a number of ways. 
The first one is to employ spectroscopic absorption feature in- 
dices, like the Lick/IDS ones (Burstein et al. |T984 Trager et 
al. 1998). Balmer line indices are strongly sensitive to the pop- 
ulation age, whereas metal line indices provide estimates of 

Send offprint requests to: Maurizio Salaris 



the GC chemical composition (see, e.g., Beasley et al. [2000; 
Strader et al. 2005 and reference therein for examples of ap- 
plication of these indices to the study of extragalactic GCs). A 
second possibility is to use suitable broad-band colour-colour 
diagrams, whereby one of the two colours is more sensitive to 
the age of the population, and the other is more affected by 
the metallicity. A widely used diagram is (V - I) - {V - K) 
(see, e.g., Hempel et al. 120031 Hempel & Kissler-Patig 120041 
Larsen, Brodie & Strader 2005 ) but there are also other pos- 
sible colour combinations, that include at least one near-IR fil- 
ter as summarized and discussed in James et al. ( 2006 ) - see 
also Rejkuba (2001 ) and Hempel et al. (2005). In addition, the 
(U - B) - (B - V) diagram can in principle also be used to 
disentangle age and metallicity effects in case of old popula- 
tions (see, e.g., Yi et al. 1200 4 for an application to the GCs 
of NGC 5128). A third possibility to break the age-metallicity 
degeneracy is to fit the observed spectral energy distribution 
(SED) of the cluster - where SED denotes in this case an en- 
semble of absolute magnitudes in a given set of broad-band 
filters - as discussed by Anders et al. (2004) and De Grijs et 
al. ( 120051 see also Fan et al. 2006). In these papers the authors 
consider broad-band filters ranging from U to H, and show how 
a SED composed of the BVIH or UBVI filters is best suited to 
break the degeneracy in the GC age regime. 

The most accepted view about the observed GC colour bi- 
modality is that it is driven mainly by metallicity (but see also 
Yoon et al. 2006 for an alternative explanation based on the 
morphology of the clusters' Horizontal Branch) with the blue 
subpopulation being the metal poorer one. The age distribu- 
tion of these two subpopulations is however still a matter of 
debate. As reviewed in detail by Brodie & Strader (2006) the 
metal poor GCs are generally considered to be old, of age com- 
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parable to the age of the universe and of the oldest Galactic 
GCs, whereas conflicting results are obtained for the metal rich 
ones. A representative case of this uncertainty is the GC sys- 
tem of NGC 4365, recently reanalyzed in detail by Larsen et 
al. (|2005 1. The authors discuss how different datasets of ab- 
sorption feature indices give contradictory results about the 
presence of a metal rich young GC component; also broad- 
band colours provide conflicting results, due to offsets between 
photometries, between theoretical integrated colours obtained 
from different authors, and possible offsets between observed 
Galactic GC colours (especially (V - K)) and the theoretical 
counterparts. 

A definitive conclusion about the existence of a young com- 
ponent in extragalactic GC systems is important to constrain 
scenarios for the formation of spheroids. Necessary steps in this 
direction involve detailed assessments of systematic uncertain- 
ties related to the age determination methods using absorption 
feature indices and integrated magnitudes and colours. Here 
we will address methods based on integrated colour-colour di- 
agrams. 

A general analysis (although not specialized to the esti- 
mate of GC ages) about uncertainties in theoretical integrated 
colours of single age/single metallicity populations (also de- 
noted as Simple Stellar Populations - SSPs) was published by 
Chariot, Worthey & Bressan ( 119961 ) who compared results ob- 
tained from different sets of stellar models and different spec- 
tral libraries. More recently Yi ( 2003 1 discussed in a systematic 
way the effect of some theoretical uncertainties on selected SSP 
synthetic colours. Lamers, Anders & de Grijs (2006 1 have in- 
vestigated the effect on the photometric evolution of unresolved 
star clusters (for the case of solar metallicity) of preferential 
loss of low mass stars from a GC, due to mass segregation and 
evaporation. This process may alter appreciably the integrated 
colours of GCs. 

In this paper we focus specifically on the (V - 1) -{V - K) 
diagram used in many studies of extragalactic GCs, like the 
NGC 4365 system that - as already mentioned - provides con- 
flicting results about the cluster age distribution. In particu- 
lar, we analyze theoretically the effect of a number of system- 
atic uncertainties, on both the inferred cumulative age distri- 
bution and the [Fe/H] cumulative distribution. This follows the 
works by, e.g., Hempel et al. (2003) and Hempel & Kissler- 
Patig (2004) that show how the use of the cumulative age dis- 
tribution provides a powerful tool to establish the presence of 
a young GC subpopulation. These authors discussed in detail 
the effect of the GC sample size, observational effects due to 
contamination by background objects, and the use of differ- 
ent population synthesis models. To push further the analysis 
of this method, we investigate systematic effects related to the 
morphology of the Horizontal Branch (HB) that cannot be reli- 
ably predicted by theory yet; the metal mixture of the stars har- 
boured by the observed GCs (that cannot be inferred from the 
colour-colour diagram only); the possible efficiency of over- 
shooting from the convective cores (when present) and statis- 
tical colour fluctuations caused by the undersampling of the 
bright and fast evolutionary phases (RGB, Asymptotic Giant 
Branch) due to the low total mass of GCs. Our analysis is done 
differentially, with respect to reference synthetic samples, and 



makes use of our own stellar population synthesis database (the 
BaSTI database, see Pietrinferni et al. 120041 12006V We want to 
establish, in particular, if the effects we examine are able to 
mimic a bimodal-age distribution, even in case of a constant 
age for the whole cluster population. We also wish to assess 
eventual biases in the cluster metallicity distribution inferred 
from the colour-colour diagram, given the importance of deter- 
mining both age and metallicity distributions of GC systems to 
probe the formation histories of the parent galaxies. 

The next section presents briefly the theoretical population 
synthesis models used in our analysis, and discusses the de- 
termination of integrated magnitudes and colours. Section 3 
describes the method of construction of template GC sam- 
ples used as reference for our differential study, and Section 4 
presents a quantitative analysis of the sources of systematic un- 
certainties listed above. A summary and conclusions follow in 
Section 5. The Appendix provides useful analytical relation- 
ships to determine the cumulative age distribution of a sample 
of GCs from their (V — K) — (V - I) colours, and a test for the 
accuracy of our theoretical colour calibrations. 

2. Stellar model database and computation of 
integrated magnitudes/colours 

The stellar population model database adopted in our anal- 
ysis is the so-called BaSTI (a Bag of Stellar Tracks and 
Isochrones) that can be accessed on the web at the official site: 
http://www.te.astro.it/BASTI/index.php A description of the 
database can be found at the website as well as in Pietrinferni 
et al. (T^UOH 12006b and Cordier et al. (2006). The stellar 
model/isochrone part of the database covers 1 1 different metal- 
licities (from extreme Pop II stars to super metal rich popula- 
tions) and two different heavy element mixtures for each metal- 
licity, i.e. scaled solar and a-enhanced typical of the Galactic 
halo population. We employed consistent bolometric correc- 
tions for both scaled-solar and a-enhanced mixtures, even 
at super-solar metallicities. The evolutionary tracks follow all 
major stages until the end of the Asymptotic Giant Branch 
(AGB) or carbon ignition, depending on the value of the stel- 
lar mass. Mass loss from the stellar surface is accounted for 
using the Reimers (1975) law and two values of the free pa- 
rameter j] (77=0.2 and 0.4, respectively). Superwinds along the 
AGB phases are also accounted for. In case of stars that de- 
velop convective cores during the central H-burning phase we 
have computed models with and without overshoot from the 
Schwarzschild boundary of the central convective regions. 

We have determined the integrated magnitudes (and 
colours) presented in this work as follows. When analytical in- 
tegrated magnitudes for an SSP with a given age t and metal- 
licity Z are needed, we simply integrate the flux in the chosen 
filters along the representative isochrone; we employ the Initial 
Mass Function (IMF) by Kroupa, Tout & Gilmore (119931 to 
determine the number of objects populating each point along 
the isochrone, and an appropriately chosen normalization con- 
stant. The upper mass limit is given by the most massive star 
still evolving along the isochrone, the lower mass limit is equal 
to 0.5 M G , that is the lower mass limit of the models in the 
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BaSTI database 1 . Objects with masses below this limit con- 
tribute only negligibly to the integrated magnitudes and colours 
we are dealing with, as we have verified by implementing the 
very low mass star models by Cassisi et al. (2000 1, that extend 
down to ~0. 1M Q . The value of the IMF normalization constant 
appropriate for the specified total mass M, is obtained by en- 

r M„ 

forcing the condition M, = J M tff(M)MdM where tff(Af) is the 
IMF, M is the stellar mass, M; and M u the lower and upper 
mass values set, respectively, to 0.1 M Q (the contribution to the 
total cluster mass of objects with masses below 0.5 M Q can- 
not be neglected) and to the highest stellar mass still evolving 
in a population of that age and Z. The resulting analytical in- 
tegrated colours are independent of the total cluster mass M,. 
Whenever integrated magnitudes and colours are needed for a 
(f, Z) combination not available in the isochrone database, they 
are determined by quadratic interpolation within our (?, Z) grid. 

The analytical computation described above is somewhat 
the standard procedure followed in the literature. However, 
this procedure is strictly valid only when the number of stars 
is formally infinite. Whereas the analytical computation im- 
plies that all points along the isochrones are smoothly popu- 
lated by a number of stars that can be equal to just a fraction 
of unity in case of fast evolutionary phases, in real clusters 
the number of objects at a point along the Colour Magnitude 
Diagram (CMD) is either zero or a multiple of unity. When 
the number of stars harboured by the observed population is 
not large enough to sample smoothly all evolutionary phases, 
statistical fluctuations of star counts will arise. This means 
that an ensemble of SSPs all with the same M t , t and Z 
will show a range of integrated magnitudes (and colours) 
due to stochastic variations of the number of objects pop- 
ulating the faster evolutionary phases, that in case of inter- 
mediate/old SSPs are the upper RGB and AGB. The mag- 
nitude of these fluctuations will be larger in the wavelength 
ranges most affected by the flux emitted by RGB and AGB 
stars, typically the near-IR and longer wavelengths. This is 
the origin of the colour statistical fluctuations we will dis- 
cuss in the next section, and that have been addressed in 
the literature by, e.g., Chiosi, Bertelli & Bressan (1988 1, 
Girardi & Bica ( fT993t . Santos & Frogel ( fT997b . Cervino & 
Valls-Gabaud 12003L Cervino & Luridiana (flUoH ). Fagiolini, 
Raimondo & DeglTnnocenti (2006 1. 

To take into account the effect of these fluctuations we have 
computed integrated magnitudes and colours with a Monte- 
Carlo formalism. For a selected pair (f, Z) we have first de- 
termined the appropriate isochrone by interpolating quadrati- 
cally (when necessary) in both age and metallicity among a 
set of BaSTI isochrones. As a next step, star masses are ran- 
domly extracted following the Rroupa et al. ( 1993) IMF, until 
a specified value of M, is reached. The position of the indi- 
vidual stars along the isochrone is determined by linear inter- 
polation between the two neighbouring tabulated mass values 
(a linear interpolation is sufficient, given the large number of 
points - 2250 - used to sample the isochrone). The fluxes of 



1 We are presently working to extend the BaSTI database to the 
regime of very low mass stars 



all synthetic stars populating the isochrone are then added, and 
magnitudes plus colours computed. 

As a technical detail, we add that to compute the value of 
the SSP mass consistently with the analytical procedure, the 
lower limit for the random extraction of stellar masses is set to 
O.1M . Whenever a mass value between this limit and 0.5M Q 
is extracted, it is added to the total mass budget but disregarded 
for the isochrone population. 

3. Reference GC populations 

As stated in the Introduction, the aim of this paper is to in- 
vestigate the effect of a number of systematic uncertainties on 
GC age and [Fe/H] cumulative distributions, as inferred from 
the integrated (V — T) — (V - K) diagram. As a first step we 
need to create some reference GC population that will be used 
for differential comparisons. We have therefore first created a 
synthetic GC system (labelled 'template' system) with the fol- 
lowing characteristics: 

- We employed a-enhanced isochrones (with a mean 
[<?/Fe]~0.4, the exact individual abundances of the metal 
distribution are listed in Pietrinferni et al. 2006 1 without 
overshooting beyond the Main Sequence convective core 
boundaries (for the appropriate mass range) and mass loss 
parameter 77=0.2 constant for all RGB stars. The HB part 
of the isochrone is a clump of stars, because of the con- 
stant value of 77. This is fairly typical of the isochrones 
used in population synthesis modelling (see, e.g., Girardi 
etal.|2000i. 

- The individual GC total masses (M,) follow the 
McLaughlin ( 1994 1 mass spectrum. This is of the 
form N(M,) oc M~p ', where we set y to 0.15 when the clus- 
ter mass is below or equal to 1.5 x 1O 5 M , and y=2.0 for 
larger masses. This is a good assumption to approximate 
the GC luminosity function (GCLF) of the Milky Way 
globulars (McLaughlin 1994 1. We considered a mass range 
of the GC population from log(M,/M Q )=3.8 up to 6.3. 

- Bimodal [Fe/H] distribution. The two components are 
Gaussians centred around [Fe/H]=-1.55 (lcr dispersion 
equal to 0.30 dex - the metal poor subpopulation) and 
[Fe/H]=-0.55 (lcr dispersion equal to 0.20 dex - the metal 
rich subpopulation) respectively. This distribution is an ap- 
proximation to the metallicity distribution of Galactic GCs 
included in the catalogue by Harris 119961) . Our template 
sample contains 240 metal poor objects and 80 metal rich 
objects. This ratio between the two populations is approxi- 
mately the same as in the Harris ( 1996) Galactic GC cata- 
logue. The total number of objects is taken as representative 
of the larger extragalactic GC samples with available mul- 
ticolour photometry. 

- Constant GC age equal to 12 Gyr. 

- Gaussian photometric errors with a lcr dispersion equal to 
0.03 mag in all filters have been applied to the individual 
magnitudes. This uncertainty is comparable to the average 
errors in the integrated colours of Galactic GCs listed by 
Harris ( 1996), and it is taken as representative of an accu- 
rate GC integrated photometry. 
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The integrated colours and magnitudes of the template popula- 
tion have been computed as follows. First, a cluster metallicity 
is extracted randomly according to the assumed input [Fe/H] 
distribution, and the age is assigned (in this case it is the same 
for all clusters). The mass M, of the cluster is then randomly 
selected according to the specified mass distribution (we do 
not introduce any correlation between the cluster mass and its 
[Fe/H]). Integrated magnitudes and colours are computed an- 
alytically as discussed in the previous section, and then per- 
turbed by a photometric error drawn randomly according to the 
specified Gaussian error law. 

The upper panel of Fig.^displays the integrated (V - I) - 
(V - K) (Johnson V- and K-, Cousins / filter) diagram of one 
single realization of the template population. Overplotted is 
also a subset of the grid of analytical integrated colours used to 
build the synthetic template. As already mentioned, our anal- 
ysis will focus on the age and metallicity distributions that 
are obtained from the position of the clusters in this colour- 
colour diagram, widely used in the literature. To this purpose, 
we stress here some relevant properties of the theoretical colour 
grid displayed in Fig.^ that will be helpful to interpret the re- 
sults presented in the rest of the paper. As a first approximation 
one can assume that (V - K) is sensitive essentially to [Fe/H], 
given that the {V-I) range spanned by lines of constant age and 
varying [Fe/H] is much smaller than the corresponding (V - K) 
range (see Fig.Q. The (V - I) colour is sensitive mainly to age 
at the lowest [Fe/H] values; however, when moving towards the 
higher end of the [Fe/H] range, the lines of constant [Fe/H] and 
varying age span also a sizable range of (V - K) colours. 

The theoretical grid covers a (V - K) interval ~2-3 times 
larger than the corresponding (V - /) range. This implies that 
similar changes of (V - K) and (V-I) colours will affect 
comparatively more the age determination (largely sensitive to 
(V-I)) than the [Fe/H] estimates. It is also important to notice 
that in the low [Fe/H] regime the lines corresponding to con- 
stant ages above 9 Gyr lie closer than at higher [Fe/H]. This 
implies that small changes of (V - I) cause larger variations of 
the retrieved age of old populations, compared to the case of 
higher [Fe/H]. 

Coming back to Fig. an obvious point is to notice the 
important effect of the photometric error, even with a lcr er- 
ror as small as 0.03 mag. Some objects appear younger than 
2 Gyr when dated according to their location in the diagram, 
and many clusters appear older than 15 Gyr. Also the individ- 
ual metallicities retrieved are obviously affected by the photo- 
metric errors, but the effect on the age is the most striking. This 
shows clearly that age-dating only a handful of objects in an ex- 
tragalactic GC population can lead to erroneous results about 
the age distribution of the cluster system. Even when larger 
samples are observed, it is paramount that their colour-colour 
distribution is not biased compared to the case for the whole 
system. 

As a second reference sample used in the comparisons that 
follow in the next section, we have created a synthetic GC sys- 
tem with an intrinsic bimodal-age population (labelled 'popu- 
lation A'). This system has the same properties of the template 
sample, the only difference being that the metal rich component 
is 4 Gyr instead of 12 Gyr old. The integrated colour-colour di- 




(V-K) 

Fig. 1. Integrated colours of clusters belonging to one single re- 
alization of the synthetic template population discussed in the 
text (upper panel). The lcr colour error bars are shown in the 
upper left corner of the panel. Also displayed is the reference 
grid of analytical integrated colours for selected metallicities 
and ages. Solid lines are constant [Fe/H] lines, dashed lines are 
constant age lines. Moving towards increasing (V - I) colours, 
the reference ages are equal to 2, 4, 6, 9, 12 and 15 Gyr. Moving 
towards increasing (V -K) colours, the [Fe/H] values are equal 
to -2.62, -2.14, -1.62, -1.31, -1.01, -0.70, -0.29, -0.09, 
0.05. The lower panel displays one single realization of the 
bimodal-age population A discussed in the text. 



agram of one realization of this bimodal-age population is dis- 
played in the lower panel of Fig.^ the analytical colour grid 
used to create this sample A is obviously the same employed 
for the template population, and it is also plotted in the dia- 
gram. To assess quantitatively the differences in the age distri- 
butions retrieved from the colour-colour diagram, we employ 
the cumulative age distributions - introduced by, e.g., Hempel 
et al. (2003 2004 1 to establish statistically the existence of 
young subpopulations - as obtained from the individual cluster 
colours and the reference grid displayed in Fig.^ The cumu- 
lative age distribution (CAD - number fraction of clusters with 
age larger than f, as a function of f,) for one single realization of 
a synthetic cluster populations has been computed as follows. 
We associate to each cluster an age greater than f, when it lies 
above the line of constant f, (we considered f,- values equal to 0, 
2, 4, 6, 9, 12 and 15 Gyr) in the (V-I)-(V-K) diagram. Objects 
below the 2 Gyr line are assigned to the bin corresponding to 
t > tj=0 Gyr. Given the (V - K) of an individual synthetic GC, 
we determine by linear interpolation within the reference grid 
the expected (V - 1) of the two neighbouring constant age lines 
at the cluster (V - K). After this, it is straightforward to use 
the cluster (V - I) to assign the object to the appropriate CAD 
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t, (Gyr) 

Fig. 2. CADs for the 12 Gyr template sample. The solid line 
displays the real age distribution, dotted and dashed lines the 
CAD retrieved from the colour-colour diagram for the full sam- 
ple (dashed line) and for only the clusters with (V - K) >2.4 
(dotted line - see text for details). 



age intervals. For the very few synthetic clusters - number of 
the order of unity - that appear in some realization beyond the 
red and blue (V - K) edges of the theoretical calibration, we 
extrapolate the constant age reference colours by fitting the last 
three points of a given f, sequence with a function of the form 
(V-I)t, - a+bln[(V-K)]. The individual CAD number counts 
are then normalized to the total number of objects in each real- 
ization. In the Appendix we provide accurate fitting formulae 
(V - 1),. = f[(V - K}] that can be applied to the whole (V - K) 
range spanned by our theoretical colours, and allow a fast but 
still reliable estimate of the age distribution from the measured 
(V - K) colours, without the need to interpolate within the ref- 
erence colour grid. 

Figures |2] and 13 show the CADs obtained from 30 realiza- 
tions of, respectively, the template and the bimodal-age sam- 
ples discussed above. Filled circles represent the mean of the 
results obtained from the multiple realizations, error bars show 
the lcr dispersion around these mean values. The size of the 
error bars associated to the individual points of the CAD is af- 
fected - keeping everything else unchanged - by the GC sam- 
ple size, and increases when the number of synthetic GCs de- 
creases (see also Hempel & Kissler-Patig 2004 1. 

Figure|2]shows very clearly the large change in the age dis- 
tribution retrieved from the colour-colour diagram, compared 
to the input one (solid line). The ages retrieved for the whole 
sample (dashed line) span a range between 2 Gyr (or less) up 
to more than 15 Gyr, and the shape of the CAD is greatly 
changed, displaying a smooth decrease with f,, because the syn- 
thetic clusters tend to be more evenly distributed over the ref- 
erence clour grid. If we restrict the sample to the more metal 
rich objects, with (V - K) > 2.4 (about half of the total sain- 




ts (Gyr) 

Fig. 3. As in Fig.|5]but for population A with the bimodal age 
distribution. The upper panel shows the CADs for the whole 
sample (real age distribution - solid line; retrieved age distri- 
bution - dashed line) whereas the lower panel shows the age 
distribution only for clusters with (V - K) >2.4 (real age distri- 
bution - solid line; retrieved age distribution - dashed line). 

pie falls into this colour interval, that corresponds to retrieved 
[Fe/H] values larger than ~ - 1.0 to -1.3, using the calibra- 
tion provided by the reference grid) the CAD is very similar, 
although not exactly identical. This is due to the different sen- 
sitivity of the colours to age (and also to metallicity) at different 
locations on the diagram. The case for the bimodal age distri- 
bution is displayed in Fig. [21 Again, the retrieved ages show a 
CAD very different from the real one; the trend with f, is, as 
expected, much smoother. Restricting the analysis the sample 
to objects redder than ( V - K) > 2.4 allows one to study a sam- 
ple with a higher ratio of old to young objects, because we are 
mainly sampling the higher metallicity and younger clusters in 
the synthetic populations 2 . 

Figure |4] displays a comparison between the retrieved 
CADs of the two populations, for the full samples (upper panel) 
and for objects with (V - K) >2.4 only (lower panel). The dif- 
ference between the two CADs is evident, in both the full and 
restricted samples. Bimodal-age CADs are always below the 
unimodal ones. This is the basis for the use of CDAs as indi- 
cators of age bimodality, as discussed in detail by Hempel & 
Kissler-Patig (2004i. The sample with (V - K) >2.4 shows a 
larger difference between the two CADs, because in this colour 
range the ratio young/old clusters in the bimodal-age popula- 

2 The lcr dispersion of the input age distribution in the age range be- 
tween 6 and 12 Gyr is due to the fluctuation of the number of objects 
populating the overlapping tails of the Gaussian metallicity distribu- 
tions, and also to the effect of the photometric error that moves objects 
belonging to the two subpopulations in and out of this colour range 
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Fig. 4. Comparison between the retrieved CADs for the 12 Gyr 
template (solid line) and the bimodal-age population A (dashed 
line). The upper panel displays the comparison for the full syn- 
thetic samples, the lower panel the comparison only for objects 
with (V - K) >2.4. 



tion is about 50%, whereas when the full sample is considered, 
it is only 25%. 

The comparison presented in Fig. 0] summarizes the main 
results of this discussion about the CAD of our two reference 
synthetic populations. Bimodal-age (young+old) populations 
display a CAD systematically lower than the case of unimodal 
(old) populations. The extent of this difference depends on both 
the age of the young subpopulation and the number ratio be- 
tween the two components. Increasing the ratio of young to 
old clusters - keeping the age distribution fixed - shifts the 
CAD towards increasingly lower values. The same happens 
when one decreases the age of the young component - keep- 
ing constant the number ratio between the two subpopulations 
- as we will show explicitely during one of the tests of the 
next section (analogous conclusions can be found in Hempel & 
Kissler-Patig |2004l l. 

Figures[5]and [6]display the cumulative [Fe/H] distribution 
(CMeD) from the same multiple realizations of the template 
and the bimodal-age A sample. As in case of the CAD, filled 
circles represent the mean of the results obtained from the var- 
ious realizations, and error bars show the l<x dispersion around 
these mean values. For each realization we associate to an in- 
dividual cluster an [Fe/H] > [Fe/H], when it lies to the right 
of the line of constant [Fe/H], predicted by the reference grid. 
We considered [Fe/H],- values equal to -2.62, -2.14, -1.62, 
-1.31, -1.01, -0.70, -0.29, -0.09, 0.05, as in Fig.[T] As for 
objects blueward of the [Fe/H]=-2.62 line - number usually of 
the order of unity - they are included into an additional bin cor- 
responding to [Fe/H]> [Fe/H], = -4.5, that is an arbitrary low 
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Fig. 5. CMeDs for the 12 Gyr template sample. The solid line 
displays the real [Fe/H] distribution, dashed lines the CMeD 
retrieved from the colour-colour diagram for the full sample 
and for only the clusters with (V - K) >2.4. Here, as in all 
following figures displaying CMeDs, we do not show the bin 
corresponding to [Fe/H],- > -4.5, that always has a normalized 
number count equal to 1 . 



value. The individual number counts of each CMeD are then 
normalized to the total number of objects in each realization. 
Similar to the case of the CAD, we determine the appropri- 
ate [Fe/H] bin of the individual objects by linear interpolation 
among the reference colour grid. For synthetic objects above 
the highest age of the reference grid, we extrapolate the con- 
stant [Fe/H] lines by approximating (V-K) as a linear function 
of (V - I), fitted to the last three points below the grid bound- 
ary, for each reference metallicity. For objects located below 
the 2 Gyr line we supplement the grid shown in Fig. ^ with 
the integrated colours of a 1 Gyr old population. If necessary, 
slightly younger ages are accounted for. 

The CMeDs plotted in Figs.|5]and ^demonstrate how the 
input and retrieved CMeDs show only marginal differences. In 
the case of both the full sample and the (V-K) > 2.4 one, input 
and retrieved metallicity distributions are very similar for both 
age distributions. The largest difference is at [Fe/H],- = -1.31 
for the full samples, where the retrieved CMeDs show counts 
larger then the intrinsic one. This is due to the fact that the 
intrinsic [Fe/H] distribution shows a dip between the peaks as- 
sociated to the two metallicity components, that falls into this 
[Fe/H] bin; the retrieved [Fe/H] values have however a much 
less pronounced dip due to the smoothing effect of the photo- 
metric error. An important result is that the different age dis- 
tribution has only a small effect on the derived CMeD. This is 
clearly illustrated by Fig.|6j where the CMeDs retrieved from 
the template sample are also displayed. 
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Fig. 6. As in Fig. [5] but for population A with a bimodal age 
distribution. For comparisons we also display - dashed-dotted 
line - the corresponding CMeDs retrieved in the case of the 
12 Gyr template. 

4. Analysis of systematic effects 

In this section we analyze to what extent systematic effects re- 
lated to the unknown morphology of the horizontal branch in 
unresolved GCs, their unknown metal distribution, colour sta- 
tistical fluctuations and the uncertain efficiency of convective 
core overshooting, are able to change the CAD retrieved us- 
ing a fixed reference grid of analytical integrated colours. As 
already mentioned in the introduction, our main goal is to es- 
tablish whether these effects are able to mimic a bimodal-age 
distribution in presence of an uniformly old stellar population. 
We also wish to assess the possible occurrence of biases in the 
cluster metallicity distribution inferred from the colour-colour 
diagram, as compared to the intrinsic one. 

The general method of analysis is the following. We create 
synthetic populations with the same f, [Fe/H], M,, photometric 
error distributions of a reference system - this latter determined 
using our reference colour grid displayed in Fig.[T] a reference 
system could be our template single-age or the bimodal-age 
population A described in the previous section - but using in- 
tegrated magnitudes that include the effect of, for example, a 
different HB morphology. Tabled summarizes the main prop- 
erties of all the synthetic GC samples employed in these tests. 

We then retrieve the age and [Fe/H] distributions of these 
new synthetic samples by employing our reference grid (that 
does not take into account, e.g., the different HB morphology) 
and compare the results with those obtained for the reference 
systems. 

The quantitative results of these comparisons are inevitably 
dependent on the assumed photometric error, [Fe/H], age and - 
in some case - M, distribution of the synthetic GCs. With the 



exception of the photometric error, that when largely increased 
tends to wash out all other effects, the general trends should 
be unaffected by the precise choices made about these param- 
eters, and provide useful guidelines to understand the impact 
of these uncertainties on our interpretation of extragalactic GC 
systems. We employed throughout our analysis the [Fe/H], t 
and M, distributions of a real GC system, e.g. the Galactic one, 
already introduced in the previous section. As for the photo- 
metric error, we continue to use a l<x=0.03 mag value, that cor- 
responds to accurate photometry of Galactic GCs. The size of 
our synthetic samples are also large - the same as the template 
and the bimodal-age ones introduced in the previous section. 
This combination of - realistic - 'small' photometric error and 
large sample size, highlights more clearly the role played by 
the uncertainties we are going to address. As mentioned pre- 
viously, the effect of reducing the sample size is purely to in- 
crease the error bars attached to the individual points of the 
retrieved CADs and CMeDs. The effect of increasing photo- 
metric errors will be addressed at the end of this section. 

4.1. Horizontal Branch morphology 

To date, the theory of stellar evolution is unable to predict from 
first principles the expected morphology of the HB in an old 
stellar population of fixed age and chemical composition. This 
arises from the lack of a theory for the stellar mass loss along 
the RGB phase. Once age and chemical composition are fixed, 
the HB morphology is completely determined by the amount 
of mass lost along the RGB phase by the progenitor objects. 
From synthetic HB modelling we know that stars along the 
RGB phase have to lose on average ~O.1O-O.2OM in order 
to reproduce the typical colour location of the HBs observed 
in Galactic GCs, and a lcr spread of the order of «O.O2M 
around this mean mass loss value has to be typically consid- 
ered to match the observed HB colour extension (see, e.g., Lee, 
Demarque & Zinn 1994 1. However, the precise values may dif- 
fer from cluster to cluster and are in principle unknown for ex- 
tragalactic unresolved GCs. 

Observational determinations of mass loss rates along the 
RGB are scarce (see, e.g. Origlia et al. |2002 1 and have not 
led yet to a mass loss theory with predicting power. The 
Reimers ( 1975 1 mass loss formula is traditionally used in stel- 
lar evolution modelling; it contains a free parameter rj, whose 
value determines the efficiency of the mass loss process (see 
also Catelan 2006 for a summary of alternative, but less used 
mass loss formulae). 

A quantitative example of the effect of different RGB mass 
loss assumptions on the HB magnitudes and colours is shown 
by Figure Here we display the magnitude- T e ff HB se- 
quence for a 12 Gyr a-enhanced isochrone and the two labelled 
[Fe/H] values (that correspond approximately to the mean val- 
ues for the two metallicity components of our synthetic GC 
populations) for 77 = 0.2 (the value adopted in our template 
GC sample) and 0.4. These two values of 77 produce a total 
RGB mass loss equal to, respectively, O.O8M and O.17M 
at [Fe/H]=-1.6; O.1OM and 0.20M o at [Fe/H]=-0.7. Larger 
mass loss implies a hotter HB, due to the smaller ratio be- 
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Table 1. Main properties of the synthetic GC samples discussed in this paper (see text for details). 



Label 


Mass + [Fe/H] 


Age 


Mass loss 


Convection 


Metal mixture 


lcr phot, error 


Colour computation 


Template 


Milky Way GCs 


12 Gyr 


77=0.2 


No-oversh. 


a-enh. 


0.03 mag 


Analytic 


A 


Milky Way GCs 


12+4 Gyr 


77=0.2 


No-oversh. 


c-enh. 


0.03 mag 


Analytic 


B 


Milky Way GCs 


12 Gyr 


77=0.4 


No-oversh. 


ff-enh. 


0.03 mag 


Analytic 


B2 


Milky Way GCs 


12 Gyr 


bimodal 77=0.2,0.4 


No-oversh. 


a-enh. 


0.03 mag 


Analytic 


C 


Milky Way GCs 


12 Gyr 


0.2 < 77 < 0.4 


No-oversh. 


a-enh. 


0.03 mag 


Analytic 


D 


Milky Way GCs 


12 Gyr 


77=0.2 


No-oversh. 


a-enh. + scaled solar 


0.03 mag 


Analytic 


E 


Milky Way GCs 


12 Gyr 


77=0.2 


No-oversh. 


a-enh. 


0.03 mag 


Monte-Carlo 


F 


Milky Way GCs 


12+4 Gyr 


77=0.2 


No-oversh. 


a-enh. 


0.03 mag 


Monte-Carlo 


G 


Milky Way GCs 


12+2 Gyr 


77=0.2 


No-oversh. 


a-enh. 


0.03 mag 


Analytic 


H 


Milky Way GCs 


12+2 Gyr 


77=0.2 


Oversh. 


a-enh. 


0.03 mag 


Analytic 


I 


Milky Way GCs 


12 Gyr 


77=0.2 


No-oversh. 


a-enh. 


0.10 mag 


Analytic 


J 


Milky Way GCs 


12 Gyr 


77=0.4 


No-oversh. 


a-enh. 


0.10 mag 


Analytic 
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Fig. 7. Comparison of the HB evolution in the magnitude- 
effective temperature plane, for 12 Gyr old isochrones com- 
puted with the mass loss parameter 77=0.2 (solid line) and 0.4 
(dashed line) respectively, and the two labelled metallicities. 
The V, I and K labels denote the photometric filter used for 
each of the three pairs of tracks represented in each of the two 
panels. 



tween total mass and He-core mass; this can be clearly seen in 
the figure, together with the fact that the effect is much larger 
at the lower metallicity. The magnitude of the Zero Age HB 
(ZAHB) location (stars spend most of their HB evolution near 
the ZAHB location) becomes fainter by a large amount in / 
(about 0.8 mag) and even more in K (about 1.5 mag) when 
77 is doubled and [Fe/H]=-1.6. As a general rule and for all 
metallicities, the size of the magnitude change depends on the 
photometric filter one considers. 

Uncertain predictions of the HB morphology in unresolved 
GCs lead to uncertainties in the predicted integrated magni- 
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Fig. 8. Comparison between the retrieved CADs for the 12 Gyr 
template (solid line) the 77 = 0.4 population B (dashed-dotted 
lines) and the bimodal 77 population B2 (dotted line). The up- 
per panel displays the comparison for the full synthetic sam- 
ples, the lower panel the comparison only for objects with 
(V - K) >2.4 (the results for population B2 are not displayed 
in this panel, because they are identical to the case of the tem- 
plate population). The corresponding CADs for the bimodal- 
age population are also displayed as dashed lines. 



tudes and colours, given that the HB phase is bright and rela- 
tively long-lived, and accounts for a non-negligible fraction of 
the cluster integrated flux. 

To give a quantitative estimate of how a wrong assumption 
about the unknown HB colour affects the ages derived from the 
(V -K)-(V - 1) plane, we first created a synthetic GC popula- 
tion labelled 'population B') with exactly the same properties 
of the 12 Gyr template one, but employing isochrones com- 
puted with ?7=0.4, i.e., with a bluer HB. We then determined the 
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Fig. 9. As in Fig.[S] but for the comparison with the 12 Gyr old, 
extended HB population C (dashed-dotted line). 

CAD for 30 realizations of this sample, using the reference grid 
of our template population (with 77=0.2). The resulting CAD is 
displayed in Fig. [8] compared to the results for the template 
population and the bimodal-age population A discussed in the 
previous section. 

When the whole cluster samples are compared, popula- 
tion B shows an age distribution very similar to the bimodal- 
age population A, whereas if we consider only objects with 
(V - K) >2.4, there is essentially no difference compared 
to the result for the template population. This is because at 
[Fe/H]< —1.6 the bluer HB morphology affects strongly the in- 
tegrated colours, that become bluer and mimic a much younger 
age if the CAD is determined from the reference (77=0. 2) colour 
grid. This effect is negligible at higher metallicities because 
the HB colour and magnitude change only slightly when dou- 
bling 77 (see Fig. |7J. Of course, if the mass loss is so high as 
to produce a much bluer HB also at [Fe/H]> —1.0, then we 
would retrieve substantially young ages also for the cluster with 
(V - K) >2.4. 

As a second experiment, we tested the influence of having 
clusters with the same [Fe/H] but different mass loss efficien- 
cies along the RGB. This effectively means that GCs with the 
same metallicity (and age) can have different HB morphologies 
(as actually observed in the Galaxy, like in case of the cluster 
pair M3-M13). To this purpose we created a 12 Gyr GC pop- 
ulation (labelled 'population B2') with the same properties of 
the template sample, but where 50% of the synthetic clusters 
with [Fe/H]< —1.1 (we saw before that changing 77 from 0.2 to 
0.4 does not have an appreciable effect on the retrieved CAD 
at higher metallicities) is originated from isochrones computed 
with 77=0.4, and the remaining 50% from isochrones with with 
77=0.2 (bimodal 77 distribution). The retrieved CAD is displayed 



in the upper panel of Fig. [3] for 30 realizations of this sam- 
ple, using the reference grid of our template population (with 
77=0.2). The CAD lies, as expected, above the case of popula- 
tion B (77=0.4 for all clusters) but it is still very close to the 
bimodal-age population A for ages above 9 Gyr. Of course, a 
reduced fraction of clusters with 77=0.4 would move the CAD 
towards the template case. 

As a third quantitative test, we have determined a synthetic 
population with an HB extended in colour (labelled 'population 
C) more similar to the HBs observed in real GCs. The mass 
distribution along the HB is assumed to have a fiat profile be- 
tween the values corresponding to 77=0.2 and 0.4. Figure|9]com- 
pares the retrieved CADs with the template and the bimodal- 
age population A. Predictably, the effect on the whole sample 
is smaller, and the CADs is halfway between the template and 
population A. Clusters with (V - K) >2.4 are unaffected. A 
very similar CAD, slightly closer to the template case, is found 
from a synthetic GC sample identical to population C, but with 
an HB morphology obtained employing a mean HB mass ob- 
tained from a 77=0.3 mass loss parameter along the RGB, and a 
lcr dispersion of O.O15M around this mean value 3 . This type 
of mass distribution produces a relationship between [Fe/H] 
and the HB morphology parameter HB, - (B - R)/(B + V + R) 
(where B, V and R are in this case the number of HB stars bluer 
than the RR Lyrae instability strip, within the strip, and redder 
than the strip, respectively. We employed the relationships by 
Di Criscienzo, Marconi & Caputo 2004 for identifying the po- 
sition of the RR Lyrae instability strip in our simulations) close 
to the mean relationship for Galactic GCs, as obtained from 
data in the Harris |1996 1 catalogue. 

The main conclusion of all these tests is that the unknown 
morphology of the HB in an unresolved GC population has 
the potential to strongly affect the cluster CADs - particularly 
at low-intermediate metallicities - and mimic an old+young 
bimodal-age population, even in the case where all clusters are 
old. 

The retrieved CMeDs for the two blue-HB synthetic pop- 
ulations B and C are displayed in Figs. 1101 andll II compared 
to the intrinsic ones. Differences with respect to the intrinsic 
CMeDs are generally small. At lower [Fe/H] below ~ -0.7 the 
retrieved CMeDs display systematically (slightly) lower num- 
ber fractions compared to the intrinsic values, and it is due to 
the effect of a larger mass loss, that affects mainly the low- 
intermediate metallicities, and causes bluer integrated colours. 
This affects not only the age estimate, but also the retrieved 
[Fe/H] values. The CMeD for the sample B2 is intermediate 
between the results for populations B and C. 

4.2. Scaled solar vs a-enhanced mixtures 

The initial metal distribution of extragalactic GCs for which 
we are able to measure only integrated colours is in principle 
unknown. One can use as a guideline the fact that the Galactic 
GC system displays an a-enhanced metal mixture (see, e.g., 



3 This additional test has been performed making use of the ex- 
tended HB model library included in the BaSTI database 
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Fig. 10. As in Fig.[5]but for the comparison with the 12 Gyr old 
population B with rj-0.4 mass loss law (see text for details). 




-3.0 -2.0 -1.0 0.0 



[Fe/H], 

Fig. 11. As in Fig.[5]but for the 12 Gyr old population C with 
the extended HB (see text for details). 

Carney 1996) but this might not be the case in other GC sys- 
tems 4 . 

We display in Fig.^lli nes of constant age in the (V - K) — 
(V - I) diagram, for both scaled solar and a-enhanced metal 

4 We notice that theoretical calibrations of colour-colour diagrams 
employed to determine GC ages are usually based on scaled solar 
isochrones 



Fig. 12. Lines of constant age in the (V-K)-(V-I) diagram for 
scaled solar (solid lines) and a-enhanced (dashed lines) metal 
mixtures, as obtained from the analytical relationships given in 
the Appendix. Ages are the same as in Fig.^and increase from 
bottom to top. 

mixtures, in the (V - K) range spanned by the a-enhanced cal- 
ibration. For ease of interpretation, the smooth constant-age 
lines have been determined using the analytical relationships 
given in the Appendix. The scaled solar colours have also been 
obtained from models in the BaSTI database and are homo- 
geneous with the a-enhanced ones. The differences in the two 
theoretical calibrations are evident. At the bluest end (i.e. at 
the lowest [Fe/H]) of the theoretical calibration the colour dif- 
ferences are small, but they increase sharply with increasing 
(V - K), e.g, [Fe/H]. The net effect is to assign too young ages 
to GCs with scaled solar metal mixtures when an a-enhanced 
calibration is used. Of course the reverse is true when determin- 
ing ages of a-enhanced GCs with a scaled solar calibration. 

To give more quantitative estimates of this effect we have 
created a synthetic sample (labelled 'population D') with the 
same properties as the 12 Gyr template one, but employing a 
scaled solar metal mixture for the metal rich component, i.e. the 
component centred around [Fe/H]=-0.55. We then retrieved 
the CADs and CMeDs of this system with our reference a- 
enhanced colour grid. Figure ^] compares the CADs for this 
whole new sample and the (V - K) >2.4 clusters only, with the 
corresponding template and bimodal-age population A ones. 
The "erroneous" metal mixture attributed to the metal rich 
component simulates a bimodal-age population, albeit with a 
smaller age difference than our 12+4 Gyr system. The effect is 
clearly enhanced when the sample is restricted to only the more 
metal rich clusters with (V - K) >2.4. 

The corresponding CMeDs are compared in Figure[2] The 
differences with respect to the intrinsic distributions are small, 
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Fig. 13. Comparison between the retrieved CADs for the 
12 Gyr template (solid line) and population D (dashed-dotted 
lines) whose metal rich subpopulation is assumed to be born 
out of a scaled solar metal distribution. The upper panel dis- 
plays the comparison for the full synthetic samples, the lower 
panel the comparison only for objects with (V - K) >2.4. The 
corresponding CADs for the bimodal-age population A are also 
displayed as dashed lines. 



almost identical to the case of the template and bimodal-age 
samples shown in Figs.[3]and|6] This means that the metal mix- 
ture does not influence appreciably the retrieved [Fe/H] distri- 
butions, at least for the two specific mixtures used in this test. 

4.3. Statistical fluctuations 

In Section 2 we already discussed the fact that when the num- 
ber of objects harboured by a stellar system is not large enough 
to sample smoothly all evolutionary phases, statistical fluctu- 
ations of its integrated colours do arise, caused by stochastic 
variations of the number of objects populating the faster evolu- 
tionary phases. 

The effect on the (V - I) - (V - K) plane is shown by 
Fig- El where we display the integrated colours obtained from 
500 realizations of a 12 Gyr old population with [Fe/H]=-1.62, 
and two different values of the total mass M t , namely M, = 
1.5 1O 5 M (that is the value corresponding to the peak of the 
GCLF) and M, = 1O 6 M . For each realization we have deter- 
mined the integrated colours using the Monte-Carlo procedure 
presented in Section 2. 

The points are centred around the corresponding analytical 
values of the integrated colours, with a large scatter. The total 
(V - K) range is of the order of =s0.6 mag for M t = 1.5 1O 5 M , 
and »0.4 mag for M t = 1O 6 M . The total (V - I) range is 
«0.2 mag for M, - 1.5 1O 5 M and «0.1 mag for M, = 1O 6 M . 
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Fig. 14. As in Fig.|5]but for the 12 Gyr old Population D with 
a scaled solar metal rich subpopulation. 
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Fig. 15. Monte-Carlo integrated colours of 500 realizations of 
a 12 Gyr old GC populations with the labeled total masses and 
[Fe/H] (see text for details). 



There is a clear correlation between the fluctuations of the two 
colours. The scatter decreases with increasing cluster mass, due 
to a better sampling of the faster evolutionary phases when the 
cluster mass (hence the total number of stars) increases. The 
(V-K) colour fluctuates more than (V - I), because the in- 
tegrated flux in the Zf-band is more affected than the /-band 
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by the AGB and bright RGB evolutionary stages, that are the 
fastest ones and more prone to statistical number fluctuations. 
Due to these fluctuations only (we did not include any photo- 
metric error in the simulations displayed in Fig.ll5t the [Fe/H] 
value retrieved for an individual cluster can be wrong by as 
much as «1 dex in the low metallicity regime. Retrieved ages 
can also be biased by several Gyr. 

Figure^]shows also the colour fluctuations associated to a 
12 Gyr GC with [Fe/H]=-0.6. It is important to notice that in 
this case the distribution of points is almost parallel to the con- 
stant age line. This means that at the upper end of the Galactic 
GC metallicity distribution, statistical colour fluctuations have 
a minor effect on the age estimates using the (V - K) - (V - 1) 
diagram. 

To study the effect of colour fluctuations on the retrieved 
age distribution of an entire GC system, instead of a single clus- 
ter - including the presence of 0.03 mag 1 cr photometric errors 
- we have determined one synthetic population corresponding 
to the template one, and one corresponding to the bimodal- 
age (12+4 Gyr) population, but this time the integrated colours 
of the individual clusters are determined following the Monte- 
Carlo procedure explained in Section 2 (these two populations 
are labelled 'population E' and 'F', respectively). The analy- 
sis of these two samples allows us to check, among others, 
whether statistical fluctuations can in principle erase the dif- 
ference between single-age and bimodal-age CADs. We stress 
that the inclusion of statistical fluctuations is in principle nec- 
essary to realistically reproduce the integrated photometry of 
observed GCs, unless they have a mass distribution shifted to- 
wards masses much larger than the Galactic counterpart. 

In Figs. ^]and^]we have compared the CADs retrieved 
for these two populations E and F, with those obtained for ex- 
actly same GC properties but with integrated colours computed 
analytically. 

The general trends are clear. When only objects with (V - 
K) >2.4 are considered, the CADs are essentially unchanged 
compared to the case of analytical colours. This is a conse- 
quence of the fact that at high metallicities the colour fluctu- 
ations do not significantly affect the individual age estimates 
(see Fig.ll5>. Very different is the case when the whole samples 
are considered. The inclusion of the more metal poor clusters 
causes a significant change of the CAD, compared to the case 
of analytical colours. Significant differences start between the 
6 Gyr and the 9 Gyr age bin. The CAD that accounts for the 
statistical fluctuations is above the CAD for the analytical case, 
with a larger percentage of objects at high ages. This behaviour 
is very similar for both the single-age (Fig. 1 1 61 and bimodal- 
age fFig. 1171 populations 

To summarize, the effect of statistical colour fluctuations on 
the retrieved CADs is negligible if only the redder clusters are 
analyzed. When the full GC samples are considered, a proper 
inclusion of statistical colour fluctuations tend to shift the re- 
trieved CADs towards higher values - for ages f, >6 Gyr - com- 
pared to the case when they are neglected. This behaviour is 
very similar for both the single-age (Fig. I16l > and bimodal-age 
(Fig. Mil populations. As a consequence, relative differences 
between single-age and bimodal-age CADs are preserved. 
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Fig. 16. Comparison between the retrieved CADs for the 
12 Gyr template (solid line) and population E, that accounts 
for statistical fluctuations of the star number counts (dashed- 
dotted lines). The upper panel displays the comparison for the 
full synthetic samples, the lower panel the comparison only for 
objects with (V - K) >2.4. The corresponding CADs for the 
bimodal-age population A with analytical integrated colours 
are also displayed as dashed lines. 



As for the [Fe/H] distribution, the retrieved CMeDs that in- 
clude the effect of colour fluctuations are displayed in Figs. 1181 
and^] compared to the intrinsic CMeDs. First, one can im- 
mediately notice that, again, the assumed age distribution does 
not play a significant role in the outcome of this comparison. 
In general, the number fractions for the full samples (with the 
exception of the point at [Fe/H]=-4.5, that by definition is al- 
ways equal to unity for both intrinsic and retrieved CMeDs) 
are typically lower than the intrinsic CMeDs when [Fe/H] is 
below ~ -1.6. This is due to a sizable fraction of metal poor 
objects with anomalously blue (V - K) colours, that are as- 
signed a [Fe/H]< -2.62. At higher metallicities the differences 
with respect to the intrinsic CMeDs tend to be smaller, and 
only slightly larger than in case of photometric errors only (see 
Figs.|5]and|6}. 

As an interesting detail, one can notice that the intrin- 
sic [Fe/H] (solid lines in Figs. and I19> of objects with 
(V - K) >2.4 can reach values even lower than [Fe/H]=-1.62. 
In fact, the intrinsic CMeD corresponding to [Fe/H], > -1.62 
displays a value below unity; if all clusters in the sample have 
an intrinsic [Fe/H] higher than this value, the corresponding 
CMeD would show a value equal to unity at this point. Why do 
we find such metal poor clusters in the red subsample? This oc- 
currence is not mainly due to the photometric error employed in 
the simulation - Figures[5]and[6]display values approximately 
equal to unity at this [Fe/H] - rather to the colour fluctuations 
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Fig. 17. Comparison between the retrieved CADs for the 
12 Gyr template (solid line) and the bimodal age population 
F, with integrated colours computed via a Monte-Carlo simula- 
tion (dashed-dotted lines). The upper panel displays the com- 
parison for the full synthetic samples, the lower panel the com- 
parison only for objects with (V - K) >2A. The corresponding 
CADs for the bimodal-age population A with computed ana- 
lytical colours are also displayed, as dashed lines. 



of the most metal poor objects, that can appear at anomalously 
blue but also anomalously high (V - K) colours (see Fig. 1151. 
As a consequence the intrinsic CMeD of the (V - K) >2.4 sub- 
samples will contain anomalously metal poor objects. Another 
consequence is that the retrieved metallicity of these anoma- 
lously red objects is much higher than their intrinsic one, when 
estimated from their position in the colour-colour plane. This 
explains why the retrieved CMeD at [Fe/H], = -1.62 displays 
a value much closer to unity, compared to the intrinsic one. 

Before closing this discussion on colour fluctuations, we 
wish to address a further important point. Observations of dis- 
tant GC systems usually can acquire sufficiently precise pho- 
tometry of only the brightest clusters, populating the bright tail 
of the GCLF Brightest clusters means the most massive ones, 
for which the effect of statistical colour fluctuations is smaller. 
To simulate this effect and check the consequences on CAD and 
CMeD, we selected from our populations E and F only objects 
with integrated My < -8.5 that corresponds to M, larger than 
~ 3 10 5 M Q . The resulting samples contain about 60 objects, 
therefore the error bars on the individual points of the retrieved 
CADs and CMeDs will be larger than in the previous case. 

Figure compares the retrieved CADs of the single-age 
population E and bimodal-age population F (that include colour 
fluctuations) restricted to My < -8.5, with the results for the 
full template and population A. The effect of statistical fluctu- 
ations on the ages is now almost negligible. 
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Fig. 18. As in Fig.|5]but for the 12 Gyr old population E that 
accounts for colour statistical fluctuations. 



\ 0.8 

CD 

£l 

£0.4 
Z 0.2 


^ 1 
3? 

\0.B 
m 
fc. 

x 

£0.4 
z 0.2 



full sample 



-3.0 




(V-K)>2.4 

bimodal + statistical fluctuations 
_i i i i I i i i i I i_ 



-2.0 -1.0 
[Fe/H], 



Fig. 19. As in Fig.[5]but for the bimodal-age population F that 
accounts for colour statistical fluctuations. 



The retrieved CMeD of the population E sample ( restricted 
to My < -8.5) is compared to the intrinsic one in Fig. (^(anal- 
ogous result is obtained for the bimodal-age sample) and shows 
at low [Fe/H] much smaller differences compared to the intrin- 
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Fig. 20. Retrieved CADs for the 12 Gyr population E (up- 
per panel) and the bimodal-age population F (lower panel), 
with integrated colours computed via a Monte-Carlo technique. 
Dashed-dotted lines represent the case of GC subsamples with 
integrated My < -8.5 (about 60 clusters for each Monte-Carlo 
realization). Solid lines display the CADs retrieved for the full 
samples corresponding to the template and population A, with 
integrated colours computed analytically (see text for details). 

sic CMeD than the case displayed in Fig. because of the 
reduced fluctuations 5 . 

As a general conclusion of this further test, colour fluctua- 
tions do not play a relevant role in the determination of the age 
and [Fe/H] distributions from the integrated (V — I) — (V — K) 
diagram, when one restricts the analysis to GCs brighter than 
M v = -8.5. 

4.4. Core overshooting 

There is an ongoing debate about the possibility that in stars 
with convective cores along the Main Sequence (i.e., with mass 
larger than ~ 1.1- 1.2M ) the motion of the convective fluid 
elements is not halted shortly beyond the boundary with the 
stable radiative region as determined by the Schwarzschild 
criterion, but it goes on for an appreciable length within the 
formally stable surrounding layers (overshooting). The extent 
of the overshooting region is the focus of many investiga- 
tions like, among the most recent ones, those by Testa et 
al. J1999I . Woo & Demarque J200H and Barmina, Girardi 
& Chiosi (2002). The presence of efficient overshooting from 
Main Sequence convective cores alters both colour and lumi- 
nosity of the isochrones Turn Off (such that the Turn Off is gen- 



5 The two distributions appear formally in agreement within the lcr 
error bars, but notice the larger errors - due to the relatively small 
number of objects - associated to the M v < -8.5 sample. 
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Fig. 21. As in the upper panel of Fig. [l9]but for the sample 
restricted to objects with My < -8.5 (see text for details). 



erally brighter and hotter at fixed age) thus affecting integrated 
colours of an SSP of fixed age and chemical composition. 

This issue is irrelevant for Galactic globular clusters, given 
that they harbour stars with mass low enough not to experience 
the onset of convective cores along the Main Sequence phase; 
however, in presence of bimodal-age populations with a com- 
ponent only a few Gyr old, this phenomenon can play a role. 

Figure [33] displays our reference a-enhanced colour grid 
shown in Fig. ^ obtained from models without overshooting, 
with the addition of the 1 Gyr old sequence. Colours predicted 
by a-enhanced models including overshooting (calculated from 
isochrones in the BaSTI database, that employ the overshooting 
treatment discussed by Pietrinferni et al. 2004 1 are the same as 
the reference grid, when the age is equal or larger than 4 Gyr. 
At younger ages convective overshooting starts to play a role. 
For the sake of comparison we display in Fig.^also the 2 Gyr 
constant age line obtained from models including overshoot- 
ing. 

Quantitative and qualitative differences between the 2 Gyr 
lines with and without overshooting are remarkable. On the 
whole the 2 Gyr overshooting colour-colour sequence lies be- 
tween the 1 Gyr and 2 Gyr non overshooting ones. More specif- 
ically, for a given age and [Fe/H] pair, the (V - /) and (V - K) 
colours of the 2 Gyr overshooting models are systematically 
bluer than the values of the reference grid at 2 Gyr; the trend 
of both colours with [Fe/H] at constant age is non-monotonic 
in the overshooting case, because of the turn to the red (in both 
(V - K) and (V - I)) occurring between the two lowest metal- 
licity points. The 1 Gyr non-overshooting line shows also signs 
of a non-monotonic behaviour at low metallicities. These are 
due to the complex interplay between the variation of Turn Off 
colours, colour extension of the blue loops during the central 
He-burning phase and brightness of the RGB phase with chang- 
ing[Fe/H]. 

To highlight further the effect of convective overshooting 
on integrated colours, Fig. 1221 shows also the colour-colour di- 
agram of two synthetic samples (80 objects each) of 2 Gyr old 
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Fig. 22. Colour-colour diagram of two 2 Gyr old synthetic GC 
populations (80 clusters each) with [Fe/H]=-0.55 ±0.20, com- 
puted using models with (black filled squares) and without 
(red filled circles) core overshooting. The reference colour grid 
of Figgis a ls° displayed, together with the 2 Gyr constant 
age line obtained from models that include overshooting (thick 
dashed line) and the 1 Gyr constant age sequence from models 
without overshooting (dotted line - see text for details). The 
[Fe/H] values decrease moving away from the point with the 
largest (V - K) value, along the constant age line. 

clusters with [Fe/H] =-0.55 + 0.20 (i.e. the [Fe/H] distribution 
of the metal rich components of all synthetic populations dis- 
cussed in this paper) and lcr photometric errors of 0.03 mag, 
computed from a-enhanced rj-0.2 models both with and with- 
out convective overshooting during the Main Sequence phase. 
The inclusion of overshooting, as expected, makes the synthetic 
clusters bluer on average in both colours. This certainly af- 
fects the retrieved CADs and CMeDs, if an inconsistent colour- 
colour calibration is used. 

For a quantitative assessment of this effect we have created 
a bimodal-age population with the same [Fe/H] distribution as 
the template sample, a 12 Gyr age for the metal poor compo- 
nent, and a 2 Gyr age for the metal rich one. Photometric er- 
rors are as in the template sample. We computed the integrated 
colours analytically using the rj-0.2 c-enhanced isochrones, 
both with (this population is labelled 'population H') and with- 
out (this population is labelled 'population G') overshooting 
for the case of the 2 Gyr clusters. 

Figure |23]displays the CADs retrieved using the reference 
non-overshooting color grid for both populations G and H, 
compared to the results for the template sample, and for the 
bimodal (12+4 Gyr) population A. 

We first discuss the comparison between the two 12+2 Gyr 
populations H and G, with and without overshooting. When 



Fig. 23. CADs retrieved using the reference non-overshooting 
color-colour calibration of Fig^foi" the 12+2 Gyr bimodal- 
age populations H (dotted line) and G (dashed line). The CAD 
retrieved for the template (thick solid line) and the bimodal- 
age population A (thin solid line without symbols) are also dis- 
played. The upper panel displays the CADs for the full sam- 
ples, the lower panel the case of objects with (V - K) >2A. 



the full samples are considered, the sample that includes over- 
shooting displays a CAD differing from the non-overshooting 
one only in the youngest bins. This is easy to understand from 
Fig.|22j the "overshooting" 2 Gyr old population is located al- 
most always below the 2 Gyr line of the reference colour grid, 
whereas the "non-overshooting" population of the same age is 
- obviously - more evenly distributed around the 2 Gyr line of 
the reference grid. This explains the appreciable difference be- 
tween the two CADs in the bins corresponding to objects with 
ages larger than 2 and 4 Gyr, respectively. Differences decrease 
with increasing age because the contribution of the young pop- 
ulation decreases sharply, and what dominates is the 12 Gyr old 
component which is unaffected by overshooting. 

When we restrict the comparison to objects redder than 
(V - K)-2.4, the differences are more extreme. In this colour 
range there are very few 2 Gyr old clusters for population H 
(see again Fig. 1221 so that the large majority of objects belong 
to the old metal poor 12 Gyr component. The CAD of popu- 
lation G shows predictably much lower number fractions, be- 
cause of the presence of 2 Gyr old objects in addition to the 
12 Gyr old component. Also, the number of clusters in this 
colour range is only about 25% of the total sample for the over- 
shooting case, whereas the same fraction is about 50% for the 
non-overshooting sample. This smaller sample size explains 
the larger error bars attached to population H CAD when com- 
pared with population G. 
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Fig. 24. CMeDs retrieved for populations G and H. Dashed 
lines and filled symbols correspond to population G, dotted 
lines and open symbols to population H. Solid lines with filled 
and open symbols are the input CMeDs for populations G and 
H, respectively. The top panel shows the comparison for the full 
samples, the bottom panel the case of objects with (V-K) >2.4. 



We now compare the CAD of these two 12+2 Gyr popula- 
tions with our 12 Gyr template sample and the bimodal popu- 
lation A with 12+4 Gyr old subpopulations. As expected, the 
CADs of both 12+2 Gyr systems have systematically lower val- 
ues at each age bin than the CADs of both the template and 
population A. In case of objects redder than (V - K)-2A the 
CAD of population H population is now very close to the tem- 
plate one, because it contains essentially only 12 Gyr old ob- 
jects. This means that an analysis restricted to this colour range 
would assign to this bimodal population an essentially constant 
and old age. On the other hand the CAD of the non overshoot- 
ing population G shows, as expected, lower number fractions 
than the template and population A, even when the analysis is 
restricted to the redder clusters. 

We can now summarize the effect of Main Sequence con- 
vective core overshooting on the retrieved CAD. When a colour 
grid that does not include core overshooting is used to retrieve 
ages of a young component where overshooting is efficient, the 
retrieved CAD is biased towards young ages in the age interval 
below ~ 4 Gyr. If the analysis is restricted to the redder clusters, 
the retrieved CAD lacks almost completely the young compo- 
nent. An interesting consequence is that when one restricts the 
analysis to objects redder than a given (V - K), there is the pos- 
sibility to overlook the presence of a young component if the 
treatment of core convection in the models is not adequate. 

Turning now to the estimate of the [Fe/H] distributions, the 
CMeDs retrieved for the two bimodal-age populations G and 
H are displayed in Fig.|24] The CMeD of population G shows 



only minor deviations from the input one, whereas in case of 
population H the most evident difference is the drop between 
[Fe/H] ( - = -1.31 and -1.01 when the full sample is considered. 
This is due to the systematically lower [Fe/H] estimated for 
the metal rich 2 Gyr component (see Fig. l22> whose retrieved 

[Fe/H] has in fact a sharp drop around [Fe/H] 1 . When the 

analysis is restricted to objects redder than (V - K)-2A, pop- 
ulation H is almost completely lacking the young metal rich 
component, located at bluer colours. Compared to the case of 
population G, both input and retrieved CMeDs display a large 
drop between [Fe/H], = -1.62 and -1.01, because they are 
sampling essentially the high-metallicity tail of the metal poor 
old subpopulation. 

4.5. Varying the 1cr photometric error 

For a given M,, [Fe/H] and age distribution of the synthetic GC 
systems, the magnitude of the photometric error can affect the 
quantitative results of our comparisons (see the discussion in 
Hempel et al. 2004 1. We briefly investigate here what happens 
if we consider photometric errors larger than our assumed value 
of 0.03 mag. We examine the cases of the template 12 Gyr 
sample, and the 12 Gyr 77 = 0.4 opopulation B, discussed in 
Section l4~TI This latter sample is the one whose retrieved CAD 
appears to be the most discrepant from the template one (see 
upper panel of Fig. [8}, approximately comparable to a bimodal 
12+4 Gyr population A with 77=0.2. 

Figure [25] displays the CAD retrieved from a population 
(labelled 'population I') analogous to the 12 Gyr template, but 
including a 1 <x photometric error equal to 0. 10 mag. Compared 
to the template CAD of Fig. |3 one can notice a steeper slope 
at young ages and a flatter one at older ages. This reflects 
the much larger dispersion of points in the colour-colour di- 
agram, with sizable numbers of clusters younger than 2 Gyr or 
older than 15 Gyr. We have plotted in the same figure also the 
CADs retrieved by applying our reference grid to a population 
(labelled 'population J') analogous to the 12 Gyr GC sample 
with ?7=0.4, but including again a l<x photometric error equal 
to 0.10 mag. Comparisons with population I show that the re- 
trieved CADs are almost identical when considering both the 
whole populations or the redder clusters only. 

This result shows that lcr =0.10 mag photometric errors 
are sufficient to erase all the systematic effects (the conclusion 
for the CMeDs is similar) discussed in this paper. Given the 
similitude between the CADs retrieved for a 12 Gyr old GC 
sample with 77=0.4 and the bimodal 12+4 Gyr population with 
77=0.2 (see Section |4~T1 Fig. 1251 implies also that bimodal age 
populations with age differences of up to 8 Gyr go undetected 
from CAD analyses, in presence of photometric errors of this 
magnitude. 

5. Summary and discussion 

In the previous section we have extensively studied the effect 
of inconsistent HB morphology, metal mixture and treatment 
of core convection between a synthetic sample of GCs (corre- 
sponding to an hypothetical observed system) and the reference 
colour grid used to retrieve their age and [Fe/H]. We also inves- 
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Fig. 25. Comparison between the retrieved CADs for popula- 
tion I and J (see text and Table U for details) that include lcr 
photometric errors equal to 0. 10 mag . The upper panel displays 
the comparison considering the full cluster samples, the lower 
panel the comparison restricted to clusters with (V - K) >2.4. 

tigated the role played by statistical colour fluctuations due to 
the low M, of GCs. For all our tests we have created synthetic 
GC samples with M, (composite power law with two differ- 
ent slopes) and [Fe/H] (bimodal metal rich + metal poor com- 
ponents) distributions appropriate for the Galactic GC system, 
and a 1 <x photometric error typical of the integrated colours of 
Galactic globulars. The photometric error only is able to trans- 
form a constant or bimodal age distribution into an age distribu- 
tion that covers smoothly a larger age range, as shown already 
by Hempel & Kissler-Patig ( 2004 1 . In addition to this smooth- 
ing, due purely to photometric errors, the results of our analysis 
can be summarized as follows: 

- All effects studied in this paper may modify, to various de- 
grees, the GC age distributions obtained from integrated 
(V — K) — (V - I) diagrams. The retrieved [Fe/H] distri- 
butions are overall less affected. 

- An HB morphology bluer than the one adopted for the ref- 
erence colour grid used to estimate the GC ages can sim- 
ulate a bimodal-age population (old + young components) 
even if all GCs are coeval. This is true both in case of a 
bluer morphology for the whole cluster sample, or when 
a substantial fraction of clusters has a bluer HB than ac- 
counted for in the theoretical colour grid. For the mass loss 
laws employed in our analysis this effect is negligible when 
only the red clusters are considered (in our case clusters 
with (V - K) >2.4). More extreme amounts of mass lost 
along the RGB may influence also the ages of red clus- 
ters. Biases in the retrieved [Fe/H] distributions compared 
to the intrinsic ones are small, and restricted mainly to low 



[Fe/H], with an anomalously high number of objects as- 
signed [Fe/H] values below ~ -2.0. 

- In the presence of a scaled-solar metal rich GC component 
(the effect of changing the metal mixture is larger with in- 
creasing [Fe/H]) the use of an a-enhanced metal mixture 
in the reference colour grid introduces an age bias, but the 
retrieved [Fe/H] distribution is unaffected. In our tests the 
retrieved age distribution tends to simulate a bimodal-age 
sample - to a lesser degree than the blue HB tests - even if 
all GCs have the same age. The effect is more pronounced 
when subsamples containing only the redder (more metal 
rich) clusters are considered. 

- The effect of statistical colour fluctuations, if unaccounted 
for in the age estimate, does not affect the age distribution 
of red clusters; it tends overall to increase the number of ob- 
jects with retrieved ages larger than 6-9 Gyr, when also the 
metal poor clusters are considered. The estimated [Fe/H] 
distributions are affected to different degrees, whether the 
whole sample or only red clusters are considered. In the 
former case the metal poor clusters tend to be assigned an 
underestimated [Fe/H], whereas in the latter case there is 
the tendency to overestimate [Fe/H]. When M, is above 
« 3 1O 5 M (integrated M v > -8.5) the effect of colour 
fluctuations on the retrieved age (and [Fe/H]) distribution 
of the synthetic GC samples is almost negligible. 

- The effect of Main Sequence core overshooting is relevant 
only for ages below ~4 Gyr. If a colour grid that does not 
include core overshooting is used to retrieve ages of an 
old+young composite population where overshooting is ef- 
ficient, the resulting age distribution of the young compo- 
nent will be biased towards low ages. Also the retrieved 
[Fe/H] values will be lower than the true ones (the opposite 
is of course true if the reference grid includes overshoot- 
ing but this is not efficient in the stars belonging to the GC 
population). Subsamples containing only the redder clus- 
ters could lack almost completely the young component. 

- If the lcr photometric error on the cluster photometry 
reaches ~0. 10 mag, all systematic effects found in our anal- 
ysis are erased. Photometric errors of this size also prevent 
the detection of bimodal (old + young) age populations 
with age differences of ~8 Gyr. 

We conclude this section with a few additional considera- 
tions on GC age distributions obtained from integrated colour- 
colour diagrams. One of our conclusions is that lcr photometric 
errors of the order of 0. 1 mag would prevent the identification 
of bimodal populations with age differences of ~8 Gyr. A way 
to mitigate this problem is to make use of colours with a wider 
dynamical range. The ({/—/) - (V — K) diagram would greatly 
improve in this respect (see Hempel & Kissler-Patig 2004}. To 
this purpose, in the Appendix we provide accurate fitting for- 
mulae (17 - I) ti = f[(V - K)] (for the 77=0.2 grid) that allow 
a fast estimate of the age distribution from this colour-colour 
diagram. 

However, the use of the (U - I) colour instead of (V - I) 
does not avoid the other systematics discussed in this paper, in 
particular the thorny problem of the unknown HB morphology 
in unresolved systems. 
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Fig. 26. Comparison of the SED (UBVRIJHK bands) of three 
synthetic GCs, determined from a-enhanced isochrones with 
the labelled age, metallicity and 77 values. The integrated mag- 
nitudes are normalized in such a way that My = -6.5 (an arbi- 
trary value) for all three cases. No random photometric errors 
are included. 

Making use of both ( U - I) - ( V - K) and ( V - I) - ( V - K) di- 
agrams might, under certain conditions, help to detect the pres- 
ence of an HB population with a different morphology com- 
pared to the calibrating colour grid. The reason is explained by 
Fig- ED Here we show the SED (integrated broadband magni- 
tudes from U to K) for three different synthetic GCs computed 
with varying ages, [Fe/H] and mass loss parameter 77. After nor- 
malization to the same (arbitrary) value of My, one can notice 
the following. First of of all, all magnitudes between B and K 
are identical for a 12 Gyr [Fe/H]=- 1 .62 population with rj-QA, 
and a 7 Gyr, [Fe/H] =-1.84 population with 77=0.2. As a conse- 
quence, all colour combinations using filters from B to K (i.e. 
the (V — T) — (Y - K) diagram) show the same degeneracy, and 
the 12 Gyr ?7=0.4 population is assigned a 7 Gyr age from a 
(V-I)-(V-K) analysis that employs the reference 77=0.2 grid. 
Moreover, when 0.03 mag lcr random photometric errors em- 
ployed in our tests are applied to the magnitudes of the 12 Gyr 
system, the (V - I) and (V - K) colours of the 4 Gyr, 77=0.2 
population (and even lower ages) are well within reach. This 
kind of comparison explains the results displayed in Fig.|8] 

The normalized U magnitude is however different between 
the 12 and the 7 Gyr populations, implying that colours in- 
volving the £/-band (i.e. (U - I)) are not coincident. The con- 
sequence for the age estimates is that if one employs both 
(U - I) - (V - K) and (V - I) - (V - K) diagrams, an in- 
consistent HB morphology between observed population and 
theoretical grid, could potentially manifest itself through dif- 
ferent age distributions obtained from the two colour-colour 
diagrams. Unfortunately, the difference between the normal- 
ized U magnitudes of the 7 Gyr and the 12 Gyr populations is 
equal to only 0.05 mag (it is of 0. 10 mag between the displayed 



4 Gyr and 12 Gyr populations). The presence of lcr photomet- 
ric errors of, e.g. 0.03 mag, tends to smear out this difference 
and would make the use of the f/-band probably less helpful. 
These figures are clearly dependent on the exact morphology 
of the HB in the unresolved populations. Larger differences - 
compared to our test - between the HB morphologies of the 
observed population and the theoretical grid may still make the 
(U-I)-(V-K)and(V-I)-(V-K) CAD analysis a very use- 
ful diagnostic of the presence of blue HB objects unaccounted 
for in the theoretical analysis. 
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Appendix A: Useful analytical fits for determining 
Cumulative Age Distributions 

The CADs presented in this paper have been computed by as- 
sociating to each individual cluster an age greater than f,, when 
it lies above the line of constant f, in the (V - I) - (V - K) 
plane. Lines of constant f, for our reference a-enhanced, non- 
overshooting isochrones with rj-0.2, are well represented - 
within 0.01 mag in (V - 7) and 0.025 mag in (U - I) - by the 
following analytical relationships (the range of validity of the 
fits is also given): 

(V - 1)2 Gyr = 0.9807 ln[(V - K)] + 0.036 

(U - 1) 2 cyr = 3.1552 ln[(V - K)] - 1.097 

2.06 < (V - K)< 2.82 
(V - 7) 4 Gyr = 0.8963 ln[(V - K)] + 0.148 
(U - 7) 4 cyr = 3.1487 ln[(V - K)] - 1.022 

2.10 < (V-K) < 2.96 
(V - 1) 6 cyr = 0.8605 ln[(V - K)] + 0.198 
(U - 7) 6 Gyr = 3.2728 ln[(V - K)] - 1.106 

2.13 <(V-K)< 3.02 
(V - 1)9 cyr = 0.8390 ln[(V - K)] + 0.233 
(U-I) 9Cyr = 3.3364 h[(V-K)] - 1.132 

2.16 < (V-K) < 3.10 
(V - 7)i2 Gyr = 0.8572 ln[(V - K)] + 0.224 
(U - I) n cyr = 3.4256 ln[(V - K)] - 1.197 

2.18 < (V-K) < 3.18 
(V-7)i 5 Gyr = 0.8781 ln[(V - K)] + 0.210 
(U - 7)i5 cyr = 3.4613 ln[(V - K)] - 1.208 

2.18 < (V-K) < 3.25 

When overshooting is included, the 2 Gyr sequences dis- 
plays a non monotonic behaviour for [Fe/H]< -2.14, in both 
(V - I) and (U - I). When [Fe/H]> -2.14, the sequences are 
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well reproduced (within 0.015 mag in (V — I) and 0.03 mag in 
(17-7)) by: 

(V - I) 2 Gyr,ov = 1-0048 ln[(V - K)] - 0.018 

(17 - 7)2 = 2.5858 In[(V - K)] - 0.61 1 

1.64 < (V - K) < 2.60 

corresponding to [Fe/H] between +0.05 and -2.14. The colours 
at [Fe/H]=-2.62 are (V - K)=2M, (V - 7)=0.66 and (U - 
7)=0.98. 

Lines of constant f, for our scaled-solar non-overshooting 
isochrones with 77=0.2 are well represented (within 0.015 mag 
in (V - I)) by: 

(V - 7) 2 G yr = 0.8356 ln[(V - K)] + 0.146 
(C7 - 7) 2 Gv) - = 3.0462 ln[(V - K)] - 0.999 

2.05 < (V - K) < 2.91 
(V - 7) 4 Gyr = 0.7958 ln[(V - K)] + 0.226 
(17 - 7) 4 Gyr = 1-2806 (V-K)- 1.341 
2.11 < (V-7Q < 3.05 
(V - 7) 6 Gyr = 0.7577 ln[(V - K)] + 0.283 
(17 - 7) 6 Gyr = 1.2893 (V- K) - 1.332 
2.15 <(V-70< 3.17 
(V - 7) 9 Gyr = 0.7617 ln[(V - K)] + 0.297 
(17 - 7) 9 Gyr = 1.3492 (V-K)- 1.455 

2.17 <(V-7f)< 3.26 

(V - 7)n c,t = 0.7734 ln[(V - K)] + 0.296 
(U - 7)i2 c„- = 1-3500 (V - r> - 1.425 

2.18 <(V-70< 3.35 

(V - 7)i5 c,t = 0.8056 ln[(V - 7Q] + 0.269 

(t/ - 7)i5 c,t = 1-3617 (V-K)- 1.438 

2.17 < (V-K) < 3.43 

The 2 Gyr constant age sequences that include overshoot- 
ing display the same monotonic behaviour as for the older ages, 
and are well approximated (within 0.02 mag in (V - I) and 
0.035 mag in (U - I)) by: 

(V - 7) 2 Gyr,ov = 0-9194 ln[(V - K)] + 0.010 
(U - 7) 2 C yr,„ v = 1-2592 (V-K)- 1.402 
1.59 <(V-K)< 2.62 

We conclude this section with some further remarks about 
our reference (V - I) - (V - K) calibration detailed above. We 
have compared the [Fe/H] -(V- 1) and [Fe/H]-(V - K) empirical 
relationships obtained by Barmby et al. (2000 1 from a sample 
of Galactic GCs, with the theoretical counterpart derived from 
our template 12 Gyr sample, but with colours computed via 
a Monte-Carlo procedure, as appropriate when the spectrum 
of GC masses is considered. We recall that in this synthetic 



sample we have used a [Fe/H] distribution typical of Galactic 
GCs, and the photometric errors are consistent with the errors 
on the colours employed by Barmby et al. (2000 1. Also the GC 
mass distribution is consistent with the Galactic GCLF An age 
of 12 Gyr is typical of Galactic GC ages, as estimated by Salaris 

& Weiss J2UU21 

Barmby et al. ( 2000 1 derive from empirical [Fe/H] values 
and integrated colours the following relationships: 

[Fe/H] = (4.22 + 0.39) (V - 1) - (5.39 + 0.35) 75 objects 

[Fe/H] = (1.40 + 0.17) (V-K) - (4.62 + 0.36) 35 objects 

We considered 30 realizations of our synthetic GC sample, and 
determined from each sample slopes and zero points of the 
[Fe/H]-(V-7) and [Fe/H]-(V-7T) relationships. The mean val- 
ues and lcr spreads around the mean provide: 

[Fe/H] = (3.87 + 0.07) (V -I) - (5.14 ± 0.08) 

[Fe/H] = (1.44 + 0.03) (V-K)- (4.81 + 0.07) 

Both relationships are in agreement with the empirical ones 
within the lcr errors. Analogue relationships obtained from the 
?7=0.4 isochrones differ by much more than the lcr errors. 
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